Silicon intercalation into the graphene-SiC interface 
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In this work we use LEEM, XPEEM and XPS to study how the excess Si at the graphene- vacuum 
interface reorders itself at high temperatures. We show that sihcon deposited at room temperature 
onto multilayer graphene films grown on the SiC(OOOl) rapidly diffuses to the graphene-SiC interface 
when heated to temperatures above 1020 °C . In a sequence of depositions, we have been able to 
intercalate ^6 ML of Si into the graphene-SiC interface. 
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I. INTRODUCTION 

The interface between epitaxial graphene and bulk SiC 
plays a dominant role in both the growth and transport 
properties of graphene on SiCi"— For reasons not under- 
stood, graphene growth on the (0001) (C-face) is much 
faster than growth on the (0001) (Si-face).^ Nonetheless, 
it is widely recognized that the diffusion of Si through 
the two different SiC-graphene interfaces is responsible 
for these different graphene growth kinetics. In addi- 
tion to playing an important role in growth, the chemical 
nature of the interface has a dramatic effect on electron 
transport in monolayer graphene. Just like in growth, 
transport on the Si-face is different than on the C- 
face. Mobilities are nearly an order of magnitude higher 
on C-face graphene compared to Si-face graphene even 
though the structural quality of the films are apparently 
identical.^ 

The ability to understand and control the SiC- 
graphene interface will be a necessary step on the road 
to graphene electronics. Indeed work has already begun 
in this direction. Previous work on the Si-face of SiC has 
shown that Ge^^ and interfacial metals^ can influence 



the charge transfer from the SiC to the graphene. Sim- 
ilarly, H2 absorption has been shown to both neutralize 
the charge transfer and cause the buffer layer to become 
isolated from the SiC interface.^ These intercalants can 
cause substantial changes to the transport properties of 
the graphene. On the Si-face, H2 absorption at the inter- 
face can increase electron mobilities by 40%^ (even higher 
mobilities are achieved if the samples are gated). 

In this work we present data showing that large 
amounts of Si can be intercalated into the graphene- 
SiC interface by annealing Si deposited at the graphene- 
vacuum interface. Intercalation proceeds through het- 
erogeneous sites^^ in the graphene film beginning at ~ 
970 °C and becomes very rapid by 1020 °C . This process 
offers potentially new avenues to modify the graphene- 
SiC interface structure for new device geometries. 



II. EXPERIMENTAL 



The substrate used in these studies was an n-doped 
(n = 2 X lO^^cm"^) 6H-SiC. The sample was grown in a 
closed RF induction furnace using the Controlled Silicon 
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Sublimation (CSS) method at 1550 °C for 20 min.^ The 
sample was transported in air before introduction into 
the UHV analysis chamber. Prior to LEEM/XPEEM 
measurements the graphene films were thermally an- 
nealed at 500 °C in UHV. LEEM/PEEM measurements 
were performed at the I311-PEEM beam line at the 
MAXlab synchrotron radiation laboratory. 

ARPES and XPS measurements were made at the 
Cassiopee beamline at the SOLEIL synchrotron in Gif 
sur Yvette. The high resolution Cassiopee beamline is 
equipped with a modified Petersen PGM monochromator 
with a resolution E / /\E ^ 70000 at 100 eV and 25000 
for lower energies. The detector is a ±15° acceptance 
Scienta R4000 detector with resolution < ImeV and 
A/c^O.OlA at Jiw = 36 eV. The total measured instru- 
ment resolution is A£^<12meV. 

The Si was deposited using direct current heating of 
a Si wafer at 1150 °C . The Si deposition rate was cali- 
brated in three ways to check for consistency: (i) using 
the Si vapor pressure and the deposition geometry, (ii) us- 
ing the measured Si 2p intensity for multiple depositions 
assuming a layer by layer growth and (iii) measuring the 
attenuation of the graphene- C Is peak intensity before 
and after Si deposition. All three ways gave a deposition 
rate of 0.7±.05 ML/10 min. 



III. RESULTS 

We have performed a series of Si deposition experi- 
ments on top of a 2 ML graphene region of the graphene 
film using XPEEM. The 2 ML region was identified using 
Hibino et al's method of LEEM refiectivity oscillations.^- 
The insert in Fig. [1] shows the Si 2p core level spectra 
from the as-grown sample. We have fit the spectra with 
a spin orbit split Si 2p peak (S2) and a second peak at 
higher binding energy (HBE) we associate with a SiC 
shifted surface state similar to the one observed for the 
clean SiC (3x3) surface. The HBE peak forms after the 
first anneal and remains after further deposit and anneal 
cycles. 

When 1.4 ML of Si is evaporated at room tempera- 
ture onto the 2 ML graphene, the Si 2p peak is larger 
and shifted 1.4eV to lower Binding Energy (BE). Anneal- 
ing to 970 °C does not significantly affect the spectrum. 
However, after a 10 sec anneal to 1020 °C , the Si 2p 
intensity shifts 0.22eV to higher BE and decreases in in- 
tensity (see Fig.d]). There are four possible explanations 
for this apparent loss of Si: (i) evaporation of ~ 1 ML of 
Si from the surface, (ii) diffusion of a similar amount of 
Si to heterogeneous defect sites where 3D Si can nucle- 
ate (these sites cannot be underneath the photon beam), 
(iii) Intercalation of Si into the graphene, or (iv) diffusion 
of Si to the graphene-SiC interface. In cases (i) and (ii) 
the intensity decrease is simply Si leaving the probe area, 
while in cases (iii) and (iv) the Si 2p intensity decrease 
is caused by the electron mean free path that attenuates 
the photoelectron signal. As we will now show. Si in- 
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FIG. 1. Si 2p spectra and fits from the 2 ML graphene region. 
Spectra are shown for clean, after 1.4 ML of Si was deposited, 
and after the Si was annealed at 1020 °C for lOsec. The insert 
shows an expanded view of the graphitized spectra before Si 
deposition. The photon energy was hiy — 133eV. Solid lines 
are the fit peaks and background. 
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FIG. 2. Si 2p intensity after 5 ML of Si has deposited 
(solid line), after the film has been annealed for 20 sec at 
1020 °C (blue dots), and after an additional lOsec anneal at 
1020 °C (red open circles). The small shift of peak following 
the second anneal may results from work function change. 



tercalation into the graphene- SiC interface accounts for 
most of the lost intensity. 

Significant Si evaporation can be immediately ruled 
out because the Si evaporation rate at 1020 °C can ac- 
count for a loss of no more than ~ 0.3 ML. More im- 
portantly, once the intensity has dropped after the first 
anneal, additional heating at the same temperature has 
no further effect on the Si 2p intensity (see Fig. [2]). This 
means that after the initial anneal, the Si distribution 
remains stable. This stability also rules out substantial 
diffusion to heterogeneous sites. Some Si does initially 
diffuses to heterogeneous sites but the fact that fur- 
ther annealing at the same temperature does not change 
the Si 2p intensity indicates that this type of material 
flow quickly becomes negligable. 

We performed three similar Si deposit and anneal cy- 
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FIG. 3. /i— LEED image at 96eV from 6H C-face graphene. 
The LEED image is from the same 2 ML region where the Si 
2p data in Figs. [T] and [2] were taken. The principle diffraction 
rods from two rotated graphene sheets, are marked. The in- 
sert shows a blow up of the superstructure unit cell. The cell 
corresponds to a (v^x v^)gR6.59° reconstruction. 



cles. After each anneal, the Si 2p peak drops from its ini- 
tial deposition value and remains stable with additional 
annealing to 1020 °C . In addition, after each deposition 
and anneal cycle, the Si 2p intensity is larger than in 
the previous cycle. In other words the post-annealed sta- 
ble Si concentration is building up within the 5 jam area 
probed by the XPEEM. This buildup is either as Si in- 
tercalation between graphene sheets or Si accumulating 
at the graphene- SiC interface. 

We point out that none of the Si is intercalated between 
graphene sheets. This can be demonstrated by using ja- 
LEED and taking advantage of the unique stacking of 
C-face graphene. C-face graphene consists of sheets that 
are commensurately rotated with non-Bernal angles (i.e 
not 60°).^^'^^ This stacking leads to graphene super cells 
that can be seen in STM.^^^^^^^^ Because there are many 
rotational pairs in a typical LEED beam diameter and 
because the amplitude of the supercell vertical displace- 
ments are small (<C O.OSnm),^^ typical LEED patterns 
from the super cells are averaged over multiple weak pat- 
ters to produce a diffuse background. 

Using /i-LEED, we have selected a sub micron region 
of the sample to observe the supercell structure before 
and after Si annealing. The /i— LEED image from the 
2-layer region before Si deposition is shown in Fig. [3l 
The image shows a pattern that consists of the diffraction 
pattern from two vertically stacked but rotated graphene 
layers, gl and ^2, plus a superstructure pattern. The two 
graphene sheets are commensurately rotated by 13.17°. 
This commensurate angle gives rise to the supercell that 
can be used to index all the superlattice spots (see insert 
in Fig. [3]). The supercell in this region corresponds to 
{V57 X >/57)gR'6.59° in graphene units. The superstruc- 
ture is not due to adsorbates since the structure persists 
up to at least 1100 °C . 

Even after annealing 7 ML of Si, there is no effect on 



the superstructure pattern. If Si intercalated between 
graphene sheets, the large change in the inter layer spac- 
ing caused by inserting Si between graphene would de- 
stroy the small corrugation superlattice. After ruling out 
the other sources of the lost Si 2p intensity, we are left 
with the conclusion that after annealing, the majority of 
the deposited Si is incorporated between the SiC inter- 
face and the first graphene layer. 

The most direct evidence that significant Si interca- 
lates into the SiC-graphene interface comes from LEED 
and XPS measurements. Figure |4] shows the LEED pat- 
tern from two parts of the same sample after more than 
7 ML of Si was deposited and subsequently annealed at 
1020 °C . The image in Fig.HJ^a) is from the covered part 
of the sample were no Si was deposited. This is a typical 
LEED pattern from multilayer graphitized SiC. It shows 
both the graphene ring and, because the sample is thin, 
the diffraction pattern from the (1x1) SiC surface. Fig- 
ure [11(b) shows a pattern from the part of the sample ex- 
posed to Si. This LEED pattern only shows the graphene 
ring. Had the Si evaporated or remained on the graphene 
surface, the diffraction pattern would be the same as in 
Fig. m^a) or uniformly attenuated. Instead only the SiC 
spots have disappeared. This demonstrates that a sig- 
nificant portion of the deposited Si has diffused to the 
SiC-graphene interface where it attenuates the diffrac- 
tion pattern from the SiC, leaving the graphene pattern 
unchanged. 
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FIG. 4. LEED patterns from; (a) a masked part of the graphi- 
tized surface where Si was not deposited, and (b) from the 
same sample but at a region where 7 ML of Si was deposited 
and annealed. Note that the weak, but visible, SiC 1x1 pat- 
tern in (a) is not visible in (b). Both images were made at 
72eV. 
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A similar effect is seen in the Cls spectra before and af- 
ter Si intercalation. Figure [5] shows that the C Is signal 
from the SiC carbon disappears from the spectra after 
intercalation. The interface Si has attenuated the SiC 
carbon signal relative to the graphene carbon C Is pho- 
toelectron intensity. 
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FIG. 5. C Is spectra before Si intercalation (open circles) and 
after intercalation (closed circles). 

There is no evidence that the Si converts any signifi- 
cant amount of graphene to SiC at the 1020 °C annealing 
temperature. Figure [6fa) show an XPEEM map using 
the Si 2p core level intensity. Because of photoelectron 
attenuation, image contrast is proportional to graphene 
thickness (lighter areas are thinner). Figure [6l^b) shows 
a cut through a 2-, 3-, 4-, and 2-layer graphene region 
of the sample before and after Si deposit and anneal- 
ing. The plot shows that after Si has moved to the SiC- 
graphene interface, neither the graphene boundaries nor 
thicknesses are changing, i.e, there is no measurable con- 
version of graphene to SiC. 



IV. DISCUSION 
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FIG. 6. (a) Clean sample XPEEM image using the SiC-Si 2p 
peak at a BE of 25.8eV with photon energy hv — 133eV. The 
labels indicate the number of graphene layers determined by 
LEEM reflectivity, (b) Intensity profile across the vertical bar 
in (a) from the clean sample (blue open circles) and after ~7 
ML of Si has been deposited and ~6.3 ML intercalated into 
the interface (solid red line). 
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In order to quantitatively determine how much of the 
deposited Si diffuses to the interface we have performed 
a detailed analysis of the XPEEM and LEEM data from 
three deposition and anneal cycles. We use a model that 
assumes that most of the Si is uniformly distributed both 
at the surface and at the SiC interface. LEEM mirror 
mode images, which are sensitive to height variations on 
the surface, show no significant height variation after an- 
nealing. Since Si 3D islands on scales larger than the 
lOnm resolution of the microscope would be visible as 
surface height fluctuation, we can rule out large scale Si- 
island formation. We also allow a fraction of the Si to 
bind to C sites on the last graphene layer. We know that 
a stable Si complex is at the surface because it oxidizes 
when exposed to air (see Fig. [7j) and can only be removed 
by heating to 1450 °C . This explains the HBE peak in 

Fig.m 



FIG. 7. The Si 2p core level spectra after the 6.3 ML Si- 
intercalated sample was exposed to air at room temperature 
and subsequently heated to 1200 °C in vacuum. 



The model used to analyze our experimental results is 
shown schematically in Fig. [H In each cycle Ni (i=l,2, 
and 3) monolayers of Si where deposited {Ni =0.7 ML, 
N2 = IA ML and A^3 = 5 ML). After each anneal, Nei 
layers are evaporated, Nia = o^iNi layers have moved to 
the SiC interface and the remaining N^) = (1 — a)Ni — N^i 
are bound at the graphene- vacuum interface. As dis- 
cussed above, the model does not include intercalation 
of Si between graphene layers. The total Si intercalated 
into the SiC interface is Na = N^a- The measured Si 

2p photo-electron current after annealing, /^y , has con- 
tributions from: (i) the Si at the graphene-SiC interface 
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that is attenuated by both the Si within the layer and the 
graphene above it, and (ii) the total Si at the graphene- 
vacuum boundary. The intensity before annealing, /, 



Si ' 



depends essentially on the amount of Si deposited be- 
cause the surface Si and the graphene effectively attenu- 
ate any signal from the Si at the SiC interface. This is 
a reasonable assumption because the electron mean free 

I 



paths in graphene and Si (A^ and A^i, respectively) are 
small at the kinetic energies used in these experiments 
(see below). 



Assuming that the Si 2p cross section is the same re- 
gardless of how it is bound, we can write the relative 



intensity change after each anneal, I'si /^si-> 
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FIG. 8. (a) Model of how an Ni monolayer Si film is redis- 
tributed after annealing. Some of the surface Si has formed 
a carbide bonded to the top graphene layer, (b) Schematic 
of the surface after 7 ML have been deposited and annealed. 
Note that the surface Si and surface carbide are treated as 
part of the total surface coverage. 



where dsi and dg are the bulk Si and graphene inter- 
layer spacing, respectively. In Eq. [Talwe assume that the 
evaporated Si in the first step, A^ie, is zero because of 
the formation of a partial monolayer stable carbide as 
discussed above. We also assume that in the third depo- 
sition [Eq. [Tc] , where the Si interface layer has become 
very thick (i.e, Nadsi ^ A^^), that we can treat the in- 
terface Si as an infinite film. 



I 

The mean free path in Si at the kinetic energies used 
in these experiments (27.7eV) is 3.3A.^^ The mean free 
path in graphene at the same low kinetic energy used in 
this experiment is less well known. We have estimated 
A^ by comparing the Si 2p intensity from a 7 ML Si 
film on graphene to the the Si 2p from SiC attenuated 
by two graphene layers. This is done by assuming the 
Si film has a uniform bulk Si density and the SiC is a 
bulk terminated surface. Assuming the photo-electron 
production cross-section is the same in both layers, we 
estimate A^ = 2.8 ± 0.1 A. For comparison the mean free 
path in SiC is 3.2A.^^ We have used both values in the 
coverage calculations. The variance in the interface Si 
coverage introduced by using these two values for A^ is 
0.1 ML. 

The bulk Si evaporation rate was used as an upper 
limit estimate for Nei for both the second and third an- 
neal cycles. This gives A^ei = = 0.3 ML. This is 
small compared to the total deposited film and therefore 
does not significantly affect the total intercalation value. 
It does however influence how much Si remains on the 
surface. 

Using the experimental integrated Si 2p intensities ra- 
tios, we can solve Eq. ma)-(c) for the three cycles. The 
results are summarized in Table HI A total of 6.3 ML 
of Si was intercalated into the graphene-SiC interface. 
This leaves 0.8 ML that has either evaporated, diffused 
to defect sites, or remains on the surface. If we use the to- 
tal estimated evaporation rate of 0.3 ML, then an upper 
limit of 0.8-0.6=0.2 ML of Si remains unaccounted for 
(presumable by diffusion to heterogeneous sites). A sec- 
ond estimate of the amount of Si remaining on the surface 
can be made by using the attenuation of the graphene C 
Is core level peak before and after intercalation. This 
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TABLE I. Fit parameters for each cycle. All coverages are in 
monolayers. 



Cycle 


Deposited Si 


a,{%) 


Intercalated 




Surface 


1 


0.7 


69 


0.5 





0.2 


2 


1.4 


71 


1.0 


0.3 


0.1 


3 


5.0 


96 


4.8 


0.3 


-0.1 


Total 


7.1 


89 


6.3 


0.6 


0.2-0.5 



^ Estimated from bulk evaporation rate. 



estimate gives 0.5 ML of Si on the graphene film. This 
is an upper estimate since part of the remaining Si can 
(and does) diffuse to heterogeneous sites. 



V. CONCLUSION 

We have shown that Si rapidly diffuses from the 
graphene- vacuum surface to the graphene- SiC interface 
at 1020 °C . Up to 6.3 ML of Si were intercalated into 
the graphene-SiC interface in this study. The reason the 
Si intercalates is not known. However, at 1020 °C , the 
SiC phase diagram predicts that any excess Si will stay 
segregated. Since the excess Si has the choice of wetting 
either the graphene at the vacuum boundary or wetting 



both the graphene and the carbon rich SiC interface, the 
lower energy choice is intercalation between the SiC and 
graphene. From an equilibrium view point, the excess 
Si should bond with the excess carbon (in the form of 
graphene) to form a new stable SiC layer. The fact that 
it does not, indicates that either the surface phase dia- 
gram allows both an excess Si and C phase or there is an 
appreciable kinetic barrier to SiC formation. 

A partial Si layer (~ 0.2 — 0.5 ML) forms at the vac- 
uum interface during the first deposition and anneal cycle 
and remains tightly bound to the top layer of graphene 
during subsequent depositions and anneal cycles. This 
layer forms an oxide when exposed to air that can only 
be removed by heating above 1400 °C . There is evidence 
that above 1400 °C , Si can be removed from the inter- 
face, but more work needs to be done before this transfer 
of Si between the vacuum and the SiC interface can be 
understood. 
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